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Abstract. 

The persistence length of a single, strongly charged, stiff polyelectrolyte chain is investigated 
theoretically. Path integral formulation is used to obtain the effective electrostatic interaction 
between the monomers. We find significant deviations from the classical Odijk, Skolnick and 
Fixman (OSF) result. An induced attraction between monomers is due to thermal fluctuations 
and correlations between bound counterions. The electrostatic persistence length is found to be 
smaller than the OSF value and indicates a possible mechanical instability (collapse) for highly 
charged polyelectrolytes with multivalent counterions. In addition, we calculate the amount 
of condensed counterions on a slightly bent polyelectrolyte. More counterions are found to be 
adsorbed as compared to the Manning condensation on a cylinder. 



Polyelectrolytes (PEs) are polymers that have ionizable groups. When dissolved in so- 
lution, they dissociate into charged polymer chains and a cloud of free, mobile counterions 
carrying opposite charges jj], ^| . Such macromolecules appear in numerous industrial ap- 
plications as well as in biological systems, introducing a new kind of biologically inspired 
electrostatics ||. The spatial conformation of a single PE chain has been at the focus of 
attention of experiments @, ||, simulations HQ and theoretical models |fl PI, |To|, |Tll, |l2|, [Til , 
showing a wide range of behavior. The delicate balance between counterion entropy and long 
range electrostatic repulsion can have opposite effects on different PEs. Broadly speaking, in 
systems containing weakly charged PEs and monovalent counterions, the electrostatic repul- 
sion dominates and makes the chains stiffer |n|,[u|. On the other hand, highly charged PEs 
with multivalent counterions experience an effective attraction fl5|| which is not well under- 
stood. This attraction leads to enhanced flexibility and, in some cases, induces collapse into a 
globular conformation jl^] . For DNA macromolecules, this is known as DNA condensation [^) . 

On a mean-field level, the Poisson-Boltzmann (PB) theory predicts only intra-chain repul- 
sion Another important prediction is the Manning condensation H|l6||, where some of 
the counterions are loosely bound to the PE chain. Models going beyond mean-field theory 
take into account correlations and thermal fluctuations |jj^|l3|| . Correlations between bound 
counterions become more significant at lower temperatures, where the ions are considered to 
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be arranged in a periodic fashion similar to a Wigner crystal of electronic systems (or a cor- 
related liquid) |Q. At high temperatures, such correlations are smeared out and become less 
important, while counterion thermal fluctuations get larger and induce an attraction (sim- 
ilar to van der Waals interactions), which competes with the usual repulsion between like 



Despite an ongoing discussion, there is still no consensus on which of the above mechanisms 
is more significant in physical situations as realized in experiments Jl7[ . In this letter, we 
propose a field theoretical approach, which takes into account both correlations and thermal 
fluctuations. This allows a consistent examination of the two contributions at intermediate 
temperatures and different charge densities of the polymer. 

In order to account for counterion condensation we employ a two-phase model ||: free 
counterions in solution are in equilibrium with a one-dimensional gas of counterions bound to 
the polymer backbone. In systems without added salt, the parameter regulating counterion 
condensation is q = zl^/a where z is the counterion valency, and e/a is the monomer linear 
charge density for monomers of size a and charge e. The Bjerrum length is defined as Ib = 
e 2 /eksT, where e is the dielectric constant and ksT is the thermal energy. We will see that q 
is the significant (temperature dependent) parameter which determines the system behavior. 
On an infinite and straight cylinder Manning condensation occurs for q > 1 , and the condensed 
ions lower the average charge density on the cylinder to an effective g e g = 1 Jig ). Below we 
will show how the condensation picture changes for a semi-flexible PE modeled by a bent 
cylinder. 

We model the polymer as a semi-flexible, worm-like chain of N monomers. The persistence 
length, Z p 3> a, is a measure of the chain flexibility and only the case of a stiff, rod-like 
polymer, l p 3> L = Na, is considered. It is generally agreed that once condensation is taken 
into account, the overall charge density on the PE is small, and the effect of free ions is to 
screen electrostatic interactions [@,|l^l. The inverse Debye-Hiickel screening length is defined 
as k _1 = [47rz(z + 1)Zbc]~ , where c is the concentration of z : 1 salt. The screening length 
is assumed to be much smaller than the polymer contour length, kT 1 <C Na (as is usually 
the case in experiments), and much larger than the monomer size, k _1 3> a (ensuring the 
validity of the continuum approach employed here). We denote the spatial conformation of 
the polymer as R(s), < s < L, and the positions along the chain of the I bound counterions 
as R(si) . . . R(sj). Up to a normalization factor, the grand-canonical partition function of 
the system is 



where the path integral is a sum over all possible spatial conformations of the chain, [i is 
the chemical potential of a one-dimensional gas of bound counterions, Hq is the Hamilto- 
nian of a neutral chain with bare persistence length Iq, and H mt is the Hamiltonian of the 
screened electrostatic interaction between all charged monomers and bound counterions. The 
monomer charges are assumed to have a uniform charge density e/a along the chain, while 
the counterions are taken as point-like charges 



charges |L3|. 



J \I=0 ' i=l J0 / 



(1) 
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where UJr) = Ibb~ kt /r is the screened electrostatic interaction in units of fceT. Because Hi nt 
of Eq. (g) contains also self-interactions, all following integrations have a lower cut-off at a 
distance of order a . Continuous concentrations are introduced in the following way 



Vr) = - [ ds£(r-R(s)) 
a Jo 
I 

>(r) = R(*)). 



(3) 



»=i 



These can be substituted into the partition function Eq. (|l]) by making use of the identity 
operator which couples the discrete and continuous concentrations. This is done using the 
path integral representation of the delta function |l^,^0|. The extra complexity of this method 
is the introduction of two new auxiliary fields, denoted ip m and ip h which couple to <fi m and 
(j) h , respectively. The partition function then reads 



Z = 



JvR(s) Y[ P<W^i[R]l exp(-# cont ), 



urn 



exp [-H id - 



= exp < / d 3 r 



ds V m (R(s)) 



^e'+ l(r) nr)+i*'f 



d 3 r d 3 r' *(r)Z*(r')t/(r - r') 



z=\ 1 ~J 

—z z 



(4) 



where the vectors $ and \& and the matrix Z are introduced to simplify notations. Next, 
£b is expanded in powers of ip h and the integrations over b and ip h can be performed. The 
effective electrostatic interaction H c g and average values of the different fields can be obtained 
by comparing the partition function of Eq. (^) with that of a neutral system (e = 0) [ Bl| . 
Thus, there is no need to integrate out the polymer degrees of freedom and {ip m }- 

The method is similar to loop expansion in field theory |19fl . The chemical potential /i is set 
so that at the straight rod conformation, the average number of bound counterions per unit 



length is n 



M 



(q-l)/z% 



(zIb) ]/z, as predicted by Manning [n6[. Since the 



bound counterion phase is in equilibrium with free counterions in solution, fj, does not depend 
on the conformation of the polymer. In our model, as the PE bends, the number of bound 
counterions is adjusted accordingly in order to maintain this equilibrium. 

Fixing the chemical potential of the bound counterion gas, rather than the average density, 
has greater resemblance to experimental systems, since the ions are only loosely bound and 
not chemically attached to the polymer. Expanding Eq. ([!]) to first order in ip h (this is in 
fact a Gaussian approximation pi]), [i turns out to be the chemical potential of an ideal 
one-dimensional gas, fi = ln(n M 7T), and the total charge density does not depend on the 
conformation of the polymer. The effective interaction Hamiltonian, i? e ff,i is found to be 
simply the screened electrostatic interaction between all charges, which are uniformly smeared 
along the polymer at a constant density, 1/zIb- 
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This is just the Hamiltonian assumed by Odijk, Skolnick and Fixman (OSF) |J[10|], and their 
expression for the persistence length can be easily reproduced 

= ^0 + ^OSF = ^0 + ,,9 9, ; (6) 

where Iq is the bare persistence length of the neutral polymer backbone and Iosf is the 
electrostatic contribution. 

Corrections to this approximation are obtained through a cumulant expansion in higher 
powers of The second order expansion provides a 3-body correction to the effective 
interaction Hamiltonian 

H cSa = -e"^|(iV - zen 2 

L /•£ pL ^ ^ e -K|R(s)-R(s')| e -K|R(s)-R(s")l 

,o Jo Jo dsds ' ds "|R( S )-R( S ')MR( S )-R( S ")l (?) 

where fi deviates from its ideal gas value and satisfies 
Ln M e^ = 

fL fL -«|R( a )-R(,')| 

i + ,(i -*>ml ^ ^' | RM _ R( ,, )r m 

Note that the Manning condensation for an infinite straight cylinder is not expected to change 
substantially beyond its mean- field value as was discussed in Ref . |22l . Following the method 
used by Odijk |sj , the electrostatic persistence length l c can now be calculated using the new 
Hamiltonian, i? c ff,i + ^cff,2- For low salt concentrations (na <C 1), the persistence length can 
be expanded in powers of 1/ In na yielding 



L = l 



OSF 



9 (2- (? )-% T ^ + 0(l/[lnM] 5 
qm(Ka) 



(9) 



Equation (j^) is our main prediction and is depicted in Fig. 1 for different counterion valencies 
z = 1,2, 3 as a function of na. At low salt concentrations (na <C 1) or high q, the persistence 
length maintains the OSF kT 2 dependence, l c ~ Zosf ~ k~ 2 . We find that the electrostatic 
persistence length depends strongly on the valency of the counterions. For monovalent counte- 
rions, l c is usually positive, indicating an effective repulsion between the monomers. However, 
its value is smaller than the one predicted by OSF. Introduction of multivalent counterions 
reduces the rigidity of the PE significantly and usually l c < 0, indicating an effective attraction 
between monomers. 

Although the model assumes a rod-like PE, we speculate that a PE mechanical instability 
can be associated with Z p = Iq + l c = 0, which represents the limit of validity of rod- like 
behavior. In Fig. 1 it can be seen that with trivalent counterions, PEs with a wide range 
of bare persistence lengths, Iq, will be in a collapsed, globule-like conformation (L < 0). 
However, we note that a detailed analysis of the chain mechanical instability requires methods 
different from the persistence length prescription used in this Letter, to correctly account for 
the collapse transition. 

Using the same method, we performed an expansion to higher orders of ip h . Equation (^|), 
valid to second order, accounts for most of the deviation from the OSF result. Third and 
fourth order terms in the expansion represent only a relatively small correction to the second- 
order. This point will be discussed in greater detail in a future publication pjj. The expansion 
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procedure employed here is asymptotic. As such, the approximation improves as higher order 
terms are taken into account up to a certain order. Beyond this order the expansion diverges. 

The leading difference between our l c and ^osf is through a prefactor which depends 
on q{— z/e/a), and which may become negative at large q. The boundary between repulsive 
(positive) and attractive (negative) electrostatic contributions to the persistence length occurs 
at q = 2, although it is not a true phase transition since the total persistence length ^ p = Iq + l c 
is still positive. Taking into account higher order corrections of Eq. ([)]) may result in a shifted 
threshold of q — 2. 

In addition, it is of interest to note that q = 2 is the boundary between fluctuation- and 
correlation-dominated regimes. The length at which the average interaction between bound 
counterions is equal to k^T is z 2 Ib- This concentration of bound counterions corresponds to 
q = 2. It is expected that for low values of the bound counterion concentration, 1 < q < 2, 
where the average electrostatic interaction between ions is smaller than ksT, thermal fluctu- 
ations will dominate over the electrostatic, zero temperature ordering. For higher concentra- 
tions, q ^> 2, models relying on the Wigner crystal picture should be more appropriate. The 
boundary between the two regimes is at q = 2. As attraction is associated with correlations, 
this comparison gives a new meaning to the q — 2 threshold. The two complementary ob- 
servations show that fluctuations by themselves are not enough to induce effective attraction. 
This argument, by no means, is a rigorous proof as it does not capture precise numerical 
prefactors. 

It is now possible to compare the persistence length of Eq. ([}]) with two previous models: 
one which takes into account thermal fluctuations Jl3| and another for correlations (correlated 
liquid) |l2| ]. The comparison with the fluctuation model is rather straightforward. It consists 
of rewriting Eq. (7) of Ref. p3| in terms of Iosf and q and expanding it in the two limits g> 1 
and q > 1. In order to compare with the correlation model (at finite re), we slightly modified 
the method used by Nguyen et al jl2j to account for finite salt concentration. Equation (13) 
of Ref. (l2) describes the interaction of a counterion with a single Wigner-Seitz cell in the 
absence of added salt. We extend this result in presence of salt by including all cells within 
a radius re . Then, the obtained expression is expanded in the same two limits. The full 
derivation will be detailed in a forthcoming, longer paper 21 1. 

For the limit 5 > 1 we find that the leading term of the negative electrostatic persis- 
tence length of our loop expansion and the correlation model is the same, l e ~ — g 2 ZoSF- In 
contrast, the fluctuation model underestimates the electrostatic persistence length, yielding 
l c — Iosf/ q 2 - 

For the other limit q = 1 + Ag, Aq <C 1 we find agreement between our model and the 
thermal fluctuation model ||l3l] : l c ~ ^osf[1 — O(Ag)] while the correlation model overestimates 
it: l e ~ ^osf[ — 1 — O(Aq)]. Our expansion qualitatively accounts for both q limits. There are, 
however, differences in the prefactors. 

Furthermore, we note that g, and not just the temperature T, is the relevant parameter 
which determines the relative contributions of correlations compared to fluctuations. The 
limits of very low or high temperatures (re — > 00 or 0, respectively) are beyond the validity 
range of our model which explicitly assumes that L _1 <C re <C a -1 . Nevertheless, because the 
leading term of Eq. (||) does not depend on Ib, the electrostatic persistence length l c neither 
vanishes nor diverges in these two limits. This shows that the loop expansion method used 
indeed takes into account both correlations and thermal fluctuations. 

Having found the partition function, other thermodynamically averaged quantities can 
also be obtained. For instance, the density of the bound counterions on a bent chain, or 
alternatively, the effective q. The effective q (g c ff) is proportional to the total charge density 
on the PE chain (monomers and average concentration of bound counterions). To first order 
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in ip , we find the Manning result: q c g = 1 for all chain conformations. To second order in 
ijj h the effective charge density on the chain is smaller 

9eff = l + (q-l)—^ -4<1. (10) 

ok z In na p^ 

where p is the average cylinder radius of curvature. On a slightly bent rod, the average 
distance between condensed countcrions is smaller and correlations become more significant 
as compared to the straight-rod case. This makes condensation of counterions more favorable 
Jl^| . As the polymer bends, more counterions are able to condense, which in turn drives 
further reduction of the persistence length. Note that the derivation of Eq. ( |io|) assumes the 
bending of the chain is only a small perturbation to the straight-rod limit (p 3> L). On a 
significantly bent chain (p ~ L), the situation is completely different, and lies beyond the 
scope of our model. In the limit of a folded chain resembling a spherical colloidal particle, 
mean- field theory predicts that counterions are unable to condense at all 0. 

In order to examine the effect of the increased condensation, we look at the asymptotic 
form of Eq. @ for q = 1 + Aq, Aq <C 1 in two cases. In the first we allow the density of 
the bound counterions to be adjusted according to the equilibrium condition with the bulk 
(this is an expansion of Eq. (||) in Aq) . In the second case we add a constraint that fixes the 
density to be according to Manning for all conformations of the polymer. Expanding in Aq 
we recalculate l c for both cases 

le = ZOSF [1 + 0{Aq 2 )] 
l^° d = hsF [1 - [1/ H K a)]Aq + 0(Aq 2 )] . (11) 

The persistence length is highly sensitive to the density of the multivalent counterions. Com- 
parison between the above two expressions shows that corrections to Manning condensation 
for bent polymer chains have a substantial influence on the persistence length. 

In conclusion, we have derived an expression for the electrostatic persistence length of a 
single, stiff, strongly charged PE, taking into account both correlations and thermal fluctua- 
tions of the bound counterions. Correlations dominate for <; > 1, whereas thermal fluctuations 
for q > 1. These two mechanisms were considered separately in Ref. |l2| and Ref. p3| , re- 
spectively. The advantage of our loop expansion method is that it takes both mechanisms 
into account. As a result, it covers a large range of counterion valencies and polymer charge 
densities and offers a possible explanation for the discrepancies between the two models. 

Our results show a considerable decrease in the polymer stiffness (via its persistence length) 
for systems with multivalent counterions, z > 2. This decrease depends on the modified 
Manning-Oosawa parameter q. The estimates obtained for the collapse of semi-flexible PEs 
in presence of multivalent counterions are reasonable and are related, at least qualitatively, to 
the phenomenon of DNA condensation. The effective Hamiltonian of Eq. (Q) has the form of 
a three-body interaction. In the current model this is the main source of attraction and can 
explain the high sensitivity of the chain rigidity to the counterion valency z. 
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Fig. 1 - The electrostatic persistence length l a as function of rea according to OSF (solid line) and 
our Eq. (B) (dashed line). Valencies are specified next to each curve. The parameters chosen are: 
a = 4A, Ib = 7 A, so that q = 1.75z. The negative Z c values for z = 2, 3 indicate a possible collapse 
transition of the PE chain. 



